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We calculate the real-time-correlation function of the Sherrington-Kirkpatrick spin-glass model 
in a transverse field. Using a careful analysis of the perturbative expansion of the functional- 
integral representation, we derive the asymptotic form of the correlation function. In contrast to 
the previous works, we find for large transverse field a power-law decay of the correlation with time 
t as t~^^^ at zero temperature and t~^ at infinite temperature. At the small field region, we also 
find a significant change in the correlation which comes from the structural difference between two 
paramagnetic phases. 
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I. INTRODUCTION 

Among disordered systems, spin glasses^ are of great interest and have been investigated over decades because of the 
simplicity of the models, the rich resulting properties, and a large amount of applications to a variety of problems such 
as neural networks and information processing,^ The Sherrington-Kirkpatrick (SK) modeiS. is one of the objectives 
which has been intensively investigated and is also known as the most successful example which yields a large number 
of notable consequences by mean-field approach. 

Quantum random magnets are also the objectives of challenge and of significance due to their relevance to the 
physics of real disordered magnetic compounds at low temperature. However there is an obstacle. Combination of 
randomness and quantum fluctuation makes the analytical approach complicated. Physically it is expected that spin- 
glass state may be disturbed by quantum fluctuations at low temperature and the system exhibits quantum phase 
transitions, which was discussed in several former worksi^i^i^ 

In this paper we single out the SK model in a transverse field among a variety of quantum random magnet models 
because this is the simplest quantum model exhibiting the spin-glass phase transition and is believed to have direct 
relationship with the physics of a real disordered compound LiHo2;Yi_j;F4. We concentrate on clarifying the role 
of the quantum fluctuations for spin-glass systems similar to former works. We use the path-integral formalism for 
dealing with the quantum fluctuation. The effect of quantum fluctuation can be incorporated as (imaginary) time- 
dependent order parameters. In Ref. 0, one of the authors proposed a method for treating such fluctuation effect, 
where the time dependence of the order parameters is integrated out for deriving a renormalized effective free energy 
expressed in terms of classically defined time-independent order parameters. This approach is sufficient for extracting 
static properties of the system, e.g., the phase diagram or time-independent quantities, which can be assessed by the 
resulting effective free energy. However in this paper we look closely at the time dependence of the order parameter to 
investigate the dynamical properties of the system. Therefore we resort to another approach to assess time-dependent 
quantities. 

We analyze the local dynamical correlation function at an arbitrary temperature T = which is defined as 
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where erf is the z component Pauli-spin operator on site i and Z = Tv exp(— ) is the partition function with a 
Hamiltonian H . In the analysis of quantum systems, the calculation of the dynamical correlation functions is one of 
the main topics and a variety of analytical techniques have been developed for 11.^!? As stated above, we consider the 
SK model in a transverse field T, 
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where the averaging of the spin coupling J^j, denoted by the square brackets in Eq.([T]), is taken with the Gaussian 
distribution. 
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With regard to the behavior of dynamical correlation function there are many preceding analyticalSiiSiiiiii^ and 
numericali^ii^ works, where Eq.([T]) with imaginary time t = —it is analyzed. Their analytical results show that 
the correlation decays asymptotically as power law of r^^ at the zero-temperature critical point. Furthermore the 
authors of Refs. [ll| andl? pointed out the relevance of decay to the same behavior of time-correlation function 
in the single-impurity Kondo modcl^^ by making use of the mapping between two models, which seems to enforce 
the validity of their result. However, the numerical calculation in Ref. [l^ indicates a somewhat smaller value for the 
power-law index and some argument is required for the discrepancy. All of the former analytical results are based on 
perturbative expansion, and in this paper we reexamine this calculation by using a careful analysis of the expansion. 
We arrive at a different conclusion for the value of the power index, which stems from how to deal with multipoint 
correlation function appearing in the form of expansion. 

The paper is organized as follows. In Sec. |TT1 we calculate the imaginary-timc-corrclation function at zero temper- 
ature by using the perturbative expansion. By referring to the diagrammatic expansion method utilized for other 
quantum systems, we calculate the asymptotic form of the correlation function and make the comparison of it with 
the known results. In addition the result is confirmed by the numerical calculation. In Sec. IIIII we move on to the 
correlation at finite temperature. We investigate it numerically and give some considerations to the result, which sug- 
gest different asymptotic behaviors of the correlation at high-temperature region from the case of zero temperature. 
Section ITVl is devoted to conclusions. 



II. CORRELATIONS AT ZERO TEMPERATURE 



A. Imaginary-time formalism 

The real-time correlation ([1]) is obtained by the analytic continuation of the imaginary-time-correlation function. 



X(r) 



iTre-('3--)V-e-^^a| 
Z 



(3) 



The time-independent static part of x{'^)i X = So dTx{T)IP, was studied in detail in Ref. as a measure of quantum 
fluctuations. In order to deal with quantum operators analytically, we use the imaginary-time path-integral repre- 
sentation. Following the standard technique we introduce replica to perform the ensemble average. The correlation 
function is expressed with the path integral of normalized vector. 



x(t) = lim 



= lim 

ra-*0 
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J I?^4'^(r)4:^(0)expff]^''dr'{*$[s(")(r')]-i/[s(")(r')]} 
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where the superscript a is the replica index running from 1 to n, and the limit n — s- is taken afterward, s is the unit 
vector on the Bloch sphere and $ is the Berry phase term. Then we put the Hamiltonian of the transverse SK model 
into the above formula and take the Gaussian average with respect to the interaction . After performing the 
Hubbard-Stratonovich transformation, we have the path integral form of one-body Hamiltonian, 



X(r) = (4i)(r)4i)(0)exp 




Tl - T2 S 




(6) 



where the angular bracket denotes the path integral over s("^(r) as 
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The Hubbard-Stratonovich field xi'''i ~ '''2) {si"\Ti)s'f\T2)) is nothing but the time-correlation function ([3]) and 
is taken to be free from the replica index. Here we neglect the time correlation between different replicas, or in other 

words time-dependent spin-glass parameter g"" (ri , T2) ~ {si"^ (ti )si" ^ (T2)) {ct ^ a') is taken to be zero, which means 
we are in paramagnetic phase. Our goal is to solve this self-consistent equation for x(t). 
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FIG. 1: The pictorial representation of the two-point correlation function Z)r('!"2 — ti) = (sz(r2)sz(ri))r with T2 > ri. 



B. Self-consistent equation in linear order 

We expand Eq.® as a power series of in order to perforin the spin integration afterward. Up to the first order, 
we find the self-consistent linear equation for x(t), 



j2 rP 

X{T)^Dr{T) + — dn dT2{Dr{T,TuT2,0)-DriT)Dr{\n~T2\)}x{Tl-T2), 
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where the two-point correlation function Drd^i — ^21) 



Ti)si^^(T2))r is given by 



Dr{\Tl-T2\) = 



„^)r-2r|Ti-T2| 



,-/3r+2r|Ti-T2| 



(8) 
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We can easily find how two s^'s contribute to this propagator by drawing pictorial spin- flip representation as in Fig.[T] 
The spin pointing along the x direction is flipped two times by the spin operator CTz while propagating from t = to 
T = p. ft is necessary to go back to the original state due to the boundary condition s(0) = s{f3)J In the same way, 
for the four-point function £'r(''"i, T2, T3, = {si^\Ti)si^\T2)s'^J'\T2)s'^^\T4))Y , the spin is flipped four times and we 
obtain 



■Cr(Tl, T2, T3, T4) = Dr{Ti -T2+T3~ T^), 



(10) 



where we assume ti > T2 > T3 > t^. As we see from Fig. [1] it is important to take notice of the time ordering of the 
spin operators. 

The integral equation ([5]) can be solved by the differential equation for %('''); 



dtxir) - 8r^9,^x(r) + i6(r* - rV^)x(r) = 0, 



(11) 



which is obtained by taking derivatives of Eq. ([8]) several times and constructing the closed- form equation in terms of 
x(t). With the boundary condition it can be solved as 



xW = ^[^r+(r) + i^r_(T)], 



where r± = ^/W±VJ . At the zero-temperature limit, this function goes to 



(12) 



(13) 
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where we assume r < [3/2 in taking the Hmit [3 oo. This result shows that, due to disorder, the energy level splits 
into two as r ^ vT^iTJ- Proceeding further to the higher order correlations, we expect energy levels split into 
many and finally form a broad distribution, as we see in the following. 

However, we must note that this interpretation is valid only when T > J. Solution (I12|) changes its behavior in 
the case of smaller transverse fields, which is attributed to the phase transition between the paramagnetic and the 
spin-glass phasesi^ When F < J, the time-independent part x — Jq d,Tx{T)/l3 is nonzero even at zero temperature 
and nonperturbative analysis is required for the zero modci^ This sharp change in behavior could be an artifact of 
the perturbative expansion. Although this change survives even if we take the higher-order correlations into account 
as we show in the following, it is not clear whether it corresponds to the phase-transition point or not. For example, 
in the random energy model^^ the corresponding phase transition is of first order and it is not possible to determine 
the transition point from perturbative calculation (see Sec. IHI Cl for detailed calculations). In the present analytical 
calculation, we assume that F is large enough such that the zero mode gives no contribution. Then the correlation 
asymptotically decays to zero in time and the perturbative expansion is justified. The behavior at small F is analyzed 
numerically in Sec. IH El 



C. Asymptotic behavior at zero temperature 

We take higher-order correlations into account in the following. In general x(t) including all of such corrections 
is not tractable, but in the case of the zero-temperature limit the evaluation turns out to be considerably simpler. 
The integral in Eq.® consists of contributions from (i) ri^2 < t, (ii) ri^2 > t, and (iii) ti < t < T2 or T2 < t < ti. 
As we found above, the first-order result (|13p is short range and decays exponentially with respect to r. Then we 
expect that the asymptotic form of xi^) is dominated by contributions from (i) and (ii) since the constraint in (iii), 
Ti T2 ^ T, is harder than that in (i) and (ii), ri ^ T2, in performing double integral, which makes the contribution 
from (iii) negligible. At the limit T = 0, the contribution from (i) is written as 

j2 /" dn f dT2[g{T -Tl)x{Tl-T2)g{T2) - g{T -T2)x{Tl-T2)g{Tl)], (14) 

Jo JO 

where g{T) — exp(— 2Ft). Similarly the contribution from (ii) is given by 

J' / dT2 / dn [g{(3 - ti)x(ti - T2)g{T2 - r) - g{f3 - T2)x{ti - T2)g{T^ - r)] . (15) 

Using the relation x(t) — x(/3 — r), we find that Eq.ljTS]) is given by Eg. p^ with the replacement r by (3 — t. 
Comparing each term in above two equations we arrive at the conclusion that the main contribution comes from the 
first term in Eq. (ll4p . We have an approximated form now, 

Xir) ^ g{T) + f dTl ( dT2g{T-Tl)x{Tl-T2)g{T2). (16) 

Jo Jo 

This equation is solved as 

Xir) = i ^e-^ir+J/2)r ^ g-2(r-,//2)rj ^ (^7) 

which is a good approximation for Eq. when F ^ J. 

Higher-order contributions can be incorporated in the same way. We find 

X = 9 + J^9*X*g + J^9*X*9*X*9^ = 7^ * 9, (18) 

1- J-9*X 

where the asterisk * represents a time integral in convolution form, 

/*5(t)= r dT'f{r-T')g{T'). (19) 



This function product satisfies commutativity f * g — g * f and can be treated similar to ordinary product. Solving 
Eg. p^ we obtain 



n=0 
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where we omitted the asterisk Using the nth convohition product of g given by 

n-l 

we obtain the result for the imaginary correlation function at zero temperature, 

x(r) = -l/i(2Jr)e-2r-, (22) 
where Ii is the modified Bessel function of order 1. This function has an asymptotic form of exponential decay, 

^ ^ \/4^(Jt)3/2 ^ ' 

at r > J, which changes to a power decay t~'^I'^ at F = J. 
The analytic continuation to the real time is obtained as 

X{t) = ^^,h(2Jt)e-^^^\ (24) 

where J\ is the Bessel function of order 1. This shows t^^l"^ decay for arbitrary F as 

,25) 

In the Fourier-transformed space, the spectrum forms a semicircle with the center = 2F as 

^(^^) ^ ^/^^IE^E^e{2J - Ic. - 2F|). (26) 

It is worth stressing here that the real-time correlation always shows a power-law asymptotic decay (1251) . while the 
imaginary-time correlation (|23p exhibits a power-law behavior only at F = J. However, as we explained in Sec. IIIBl 
we must be careful about identifying the point F = J as a critical one. Actually, it is known by preceding work o^i^i^^ 
that the value of F at the critical point should be a larger one F ~ 1.5J. In order to observe the power-law decay 
numerically in imaginary-time framework, we must identify the critical point, which turns out to be a troublesome 
task. ^3 On the other hand, in the real-time framework the power-law behavior can be observed not only in the critical 
point but also in a broad range of large F where our perturbative calculation is justified. Therefore, we can check the 
universal power- law index of 3/2 in the paramagnetic phase at zero temperature by observing the semicircle form of 
the spectrum. 



D. Comparison with previous results 

Our result shows that the correlation decays in power as t"'^/^. On the other hand, the authors of Ref. have the 
conclusion of t~'^ decay via a rather general argument. Their argument is based on a perturbative expansion as we 
did in the present paper. Actually we obtained the same form (|20p as theirs, though there is a minor but crucial 
difference. The difference originates from the evaluation of the n-point correlations with n > A. For example, let us 
recall that we have the four-point correlation Dr{Ti,T2, t^jT^) = (si^-* (Ti)si^'' (r2)si^'' (T3)si^-' (''4))r obtained in Eg. pTI)) . 
In the previous works r^*^"^ the Wick theorem, 

Dr{Ti,T2,T3,T4) ^ Dr{\Ti - T2\)Dr{\T3 - T4|) + Driln - T3\)Dri\T2 - T4I) + DrUn - Ti\)Dr{\T2 - T3I), (27) 

was applied to the evaluation of this four-point correlation. However, this expression is not compatible with the exact 
one pO|) . It is evident that the integration of the spin variable is not a Gaussian one, and accordingly Eq. ((27|) is not 
justified. In order to show that the use of the Wick theorem leads to a wrong result, we examine the self-consistent 
equation in linear order ([5]). In our approach we have the solution as Eq. (|12p . and we should recall that the form 
of four-point correlation in Eq. lTUl) was used there. On the other hand, let us suppose that factorization (P7)l by 
two-point correlation is correct. Then we can perform the Fourier transformation of xi''') to find 
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FIG. 2: The correlation function x{^) at T = for iV 



w/J 

10. The solid line represents Eq. 



with r/J = 3. 



where n is an integer and 



tanh/3r (/3r)2 



pr (/3r)2 + (7rn)2 ■ 

Performing the Fourier transformation again for going back to the imaginary-time representation, we obtain 



Xir) = I 



rtanh/3r , , rtanh/?r , / 



r+ tanh /3T^ 



r_tanh/3r_ 



(29) 



(30) 



where f± — \/T^ ± TJ tanh /ST. This does not coincide with the exact result (fT2|) . which indicates that factorization 
(P7)) is not vahd. In addition, Eg. ([50]) does not satisfy the normalization condition xi''' = 0) = 1. The authors of 
Refs. 9T(J|ll|Iil introduced a La grange multiplier to cure this defect. However, we can show that the correct result is 
not reproduced even if the Lagrange multiplier is introduced. Consequently we conclude that the factorization is not 
a correct procedure in evaluating the time-correlation function. 

Suppose again that we can apply the Wick theorem to the analysis here, then the integration range of the convolution 
in Ea. (|19p is bounded by /?, not by r, which means that the time ordering is not respected. Hence we may consider 
Fourier transformation of Eq. (PO)) as 



1- Jl-4l3^J^Dl{u 



(31) 



which was found in Ref. ll^. This has the same form as Eq. ((20)l . but the time ordering is not treated properly and a 
different result from Eq. (pS|) is obtained. 

The imaginary-time correlation was investigated numerically in Ref. Ts' using the Monte Carlo method. The 
asymptotic form was obtained as t~" with a ~ 1.2, which was considerably smaller than the previously predicted 
value a = 2. Although the authors of Ref. [l^ interpreted this value as the one which is close to a = 1, we consider 
this supports our result a = 1.5, rather than the ones by former works. 
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FIG. 3: x(i) at T = and V/J — 3 obtained from the Fourier transformation of the data in Fig. [5] The points denote numerical 
data and the lines are analytical results. The solid lines are envelope curves iiz-K~^^^{Jt)~^^'^ . 



E. Numerical calculation 



In order to check the analytical result at T = we carry out numerical calculation. The method of numerical 
analysis is based on the spectral decomposition of the correlation function x(^) = / dixit) exp(ia;i) as 



E„„. e-^'^^-S (lo {E„, E„)) \{n\a!\m)\' 



(32) 



where n and m denote energy eigenstates. This expression means that components in the summations contribute to 
x(ti;) only when lo is equal to a difference of energies between two levels. At T = 0, it is clear that only non-negative 
excitations are allowed and x(a;) forms a spectrum only in the range of positive uj as we showed above. 

The quantity x(a;) is evaluated numerically by the diagonalization of the Hamiltonian. The number of sites is taken 
to be iV = 10, and the averaging is performed over more than 10000 samples. The zero-temperature results are shown 
in Fig. [2] and are compared with the analytical result pS)) . When the transverse field F is not so small, we find a 
good agreement despite the fact that Ea. (|26|) is the Fourier transformation of the asymptotic form ((24)) . We also note 
that a similar form of xi'^) has already been obtained in larger system sizes in Ref. 1 14, which supports our analytical 
result as well. 

We also show the correlation in real-time space in Fig. [3] by using numerical Fourier transformation of x(^)- The 
power index is estimated by linear fitting of In |Rex(i)| and In |Imx(i)| peaks versus Int. Due to the finite-size effect, 
the plot in Fig. [5] shows small oscillations and smooth tails, which leads to large deviations at large t in Fig. [3] For 
this reason, the fitting is carried out within the region t < 5. The power index is roughly estimated as 1.6, which is 
consistent with the analytical result, 1.5. We expect the small deviation from the analytical one also arises from the 
finite-size effect in small t region, which is relatively smaller than in large t region. 

As in Fig. [21 when the transverse field F is decreased the spectral band around a; = 2r approaches the origin. 
We also observe a small sharp peak in the band around a certain value of lu with a; < 2F. This peak grows up and 
approaches the origin with increasing system size N. We expect this peak finally reaches the origin in infinite system 
size limit. Since the spin-glass order parameter q is equal to the static part of xl^) at T = 0,— x{^ = 0) > implies 
the spin-glass phase. Therefore we can interpret the emergence of the peak as a sign of the phase transition from 
the paramagnetic phase to the spin-glass phase^i^^ Although our system size A*" = 10 is not large enough to find the 
transition point, the transition seems to happen at a value around F ~ 1.5J, which is consistent with the previous 
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FIG. 4: x{^) at T = 2.0J for = 10. 



analyses. Since the main aim of the present paper is not to determine the precise value of the transition point, we 
defer the detailed analysis to a future work. 

III. CORRELATIONS AT FINITE TEMPERATURE 
A. Numerical results 

Next we move on to the finite temperature case. When the temperature is nonzero, the negative excitations are 
allowed and we expect double peaks around uj — ±2T. These peaks finally become symmetric in the limit of T ^ oo, 
which can be understood from the spectral decomposition (|32p . The behavior at finite temperature is difficult to 
analyze since many energy levels contribute to the correlation function. In the imaginary-time path-integral formalism, 
the length of the integral path shrinks as the temperature is increased and the assumption of asymptotic decay in 
T becomes invalid. For these reasons we do not apply imaginary-time path-integral formalism to finite temperature 
case. 

We first observe the shape of x(a;) by the numerical method instead. The results at T = 2.0J and T = oo are 
shown in Figs. 3] and O respectively. As we expect, the double peaks are formed at finite temperature. When F is 
not so small, the result at T = oo can be well fitted by the function 



(33) 



with a parameter J. The result of fitting is also shown in Fig. O The parameter is taken as J = 0.76 J, irrespective 
of the value of F. This form of fitting function is interpreted in the following. 



B. Behavior at infinite temperature 

Let us consider a perturbative analysis of the spectral representation at T = oo 



(34) 
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FIG. 5: x{^) at r = oo for A'^ = 10. The result is symmetric in a; as x{~^) ~ xi^)- The solid line is the result of the fitting 
using the function in Ea. (|33|) with J — 0.76 J. 



In the limit of J = 0, the eigenstates of the unperturbed Hamiltonian Ho = — T^^-^crf are expressed by the 
eigenstates of af . When the numbers of the spins pointing in the negative and the positive x directions are k and 
N — k, respectively, the unperturbed eigenenergy is given by E^'' = —{N — 2k)T and the degree of degeneracy is 
AfCfc = Nl/kl{N — k)l. The operator erf flips the zth site spin and changes the number fc to A; ± 1, which determines 
the selection rule between the state n and m in Eq. ([Ml) . In the first order of perturbation theory, the Hamiltonian 
is diagonalized in a subspace characterized by the quantum number k. With this quantum number we obtain the 
approximate expression, 



N — l wCfc NCk + i 
k=0 /J=l 1^=1 



X i(fc + i,i/|fffifc,M)r 



(35) 



where /i and v denote the indices of degeneracy and the ranges depend on k. -EI^^ are eigenvalues at the first order 



Another form of E^j^ as a function of 



and may be expressed by linear combinations of as E^j^j^ — ^ 
Jij will not change the conclusion here. The averaging of Eq. (|35|) with respect to Jij yields a Gaussian factor as 



2r 



E^ 



lij Jij 



cx exp ■ 



N{uj - 2r)2 



(36) 



The contribution of Eq. (j35[) mainly comes from the sector with k ^ N/2, which is understood in conjunction with 
the factor 1/2^ in Eg. ([55]) . For example, in the case of fc = the number of the elements in the fii/ sum is equal to N 
and is negligible at TV ^ oo due to the multiplication by the factor 1/2^, whereas it is not suppressed when k ~ N/2. 

Although it is hard to perform explicitly the perturbative calculation even at the first order, we easily see that the 
result is expressed by a linear combination of Gaussian functions as 



X(cj) = 271 / daP{<7) 



cxp 



{uj - 2r)2 

2^2— 



exp 



{uj + 2r)2 

2(72 



(37) 



where P{(y) is the distribution function of the variance a. From this expression, we find that result (|33p can be 
reproduced by assuming the Gaussian distribution P{a) cx exp(— cr2/2 J^). The Fourier transformation of Eq. 
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nves 



cos(2rt) , , 

Xit) = (38) 

This shows that the correlation asymptotically decays as t^^. We must note that this index is obtained by a different 
mechanism as Ref. |^ where the same value of the index is concluded at zero temperature. 



C. Classical and quantum paramagnetic phase 

At finite temperature, the spin-glass phase is suppressed in increasing T and disappears at T = Jj^ The second- 
order phase-transition curve, namely, the spin-glass phase boundary, runs from T = 0, F ~ 1.5 J to T = J, F = on 
T-F plane. On the other hand, in Sec.|n]we saw the change in the behavior of at zero temperature when F is varied, 
and this change is expected to be persistent up to higher-temperature region than the spin-glass phase boundary. We 
can attribute this change to the one in moving across the boundary between the quantum paramagnetic (QP) and the 
classical paramagnetic (CP) phases, rather than between the spin-glass phase and one of the paramagnetic phases. 

The existence of the phase transition between the CP and the QP phases has already been found in the p-body 
interaction spin-glass model in a transverse field defined by the Hamiltonian 

^,......,«---<-r^cTf. (39) 

*l<*2<---<»p i 

The interaction Ji-^i^-'-ij, is random and the average is taken with Gaussian probability distribution. At p ^ oo this 
model is known as the random energy model and exhibits first-order phase transition between the CP and the QP 
phases. The imaginary-timc-corrclation function at the QP phase can be calculated from the self-consistent equation 
in terms of x(r), 

X(r) = ^4')(r)4i)(0)exp|^|]^''dri^''dr24"nn)x^-'(n -r2)4"HT2)|^ • (40) 

The factor x^~^ goes to zero at p — s- cx) unless x — ^- From this, it is clear the zeroth-order calculation becomes exact, 

/3r-2rr I -/3r+2rr 

Xir) = Dr{r) = ^^r^^ " (^1) 

The phase-transition point cannot be determined from this equation for the QP phase but can be determined by the 
static part of the correlation function for the CP phase. Considering the condition of the vanishing static part, we 
obtain a discontinuous change from x{'^) ~ 1 to Eq. ([lT|) . which also tells us the boundary between the CP and the 
QP phases reaches T — oo. When p takes a large but finite value, the phase-transition line is terminated at a certain 
finite temperature value and the transition between the CP and the QP phases turns to crossover lil 

Our numerical analysis for the SK (p = 2) model shows that the semicircle form changes to a broader one at finite 
temperature. The long tail of the spectrum reaches the origin u = Q, just as Eq. (f33|) . We also observe that a sharp 
peak around the origin, which was considered a precursor of the spin-glass transition for finite systems, is suppressed 
at finite temperature. These observations imply that the spectrum change at finite temperature indicates a crossover 
between the CP and the QP phases. This picture is supported by the analysis in Ref. where a smooth change in the 
static part of the correlation function was obtained at T > 0. However, it is known in a Langevin dynamics model 
of p-body interaction spin glass^* that the behavior at p > 2 should be distinguished from the one with p = 2. Our 
numerical calculation is carried out for rather small systems, and it is fair to say that our analysis is not conclusive 
and further studies are required to determine the properties of the finite temperature phase diagram. 



IV. CONCLUSIONS 



In conclusion, we have calculated the dynamical correlation function ([T]) in the transverse SK model Our main 
results are Eas. ([24|) and ((26l) at T = 0, and Eas. ((33|) and ([38]) at T = oo. For large transverse field, the correlation 
function asymptotically decays in time t as t~^/'^ at T = and t^^ a.tT — oo. The results of our analysis are different 
from those of previous works at T = Q^ iio, 11,12 -^.j^g pj.gggjj|; gpj^ system, it is the crucial point in the analysis that 
the simple Wick theorem cannot be applied to the multicorrelations of the spin operators. It was shown in Ref. [l^ 
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that the factorization for Gaussian variables is modified in the case of spin operators and a sign factor is introduced 
to respect time ordering of operators. We showed that the correct procedure leads to a different form of real-time 
correlation from the known one. 

In the present paper, we have analyzed the infinite range model to focus on the time dependence of the correlation 
function. It is important to study how the present result is changed at finite-dimensional models with finite range 
interactions. Spatial correlation can be related to time correlation; therefore the result of such study will be useful. 
With regard to critical exponent, there is an argument that the dynamic critical exponent z is related to the power 
index of the time-correlation function discussed in the present paper. 1^ In addition, in low-dimensional systems it 
is inferred from the general considerations^ that correlation function is significantly affected by Grifhths-McCoy 
singularities. 

Application to other random quantum systems is also an interesting future work. For the random quantum Heisen- 
berg model, the correlation function in Fourier-transformed space, xC'^)) was numerically analyzed in Refs. [13 and 
[21I . It is interesting to analyze the correlation in real-time space by using the analytical method developed here, and 
this will be our future work. 



ACKNOWLED GMENTS 



We are grateful to H. Nishimori and T. Obuchi for useful discussions and comments. One of the authors (K. Takeda) 
is supported by Grand-in-aid from MEXT/JSPS, Japan (Grant No.l8079006). 



^ M. Mezard, G. Parisi, and M.A. Virasoro, Spin Glass Theory and Beyond (World Scientific, Singapore, 1987). 
^ H. Nishimori, Statistical Physics of Spin Glasses and Information Processing: An Introduction (Oxford University Press, 
Oxford, 2001). 

^ D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 1792 (1975). 

* A.J. Bray and M.A. Moore, J. Phys. C 13, L655 (1980). 

^ B.S. Chakrabarti, A. Dutta, and P. Sen, Quantum Ising Phases and Transitions in Transverse Ising Models (Springer, Berlin, 
1996). 

® S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, 1999). 
K. Takahashi, Phys. Rev. B 76, 184422 (2007). 

* A. Georges, G. Kotliar, W. Krauth, and M.J. Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 
^ J. Miller and D.A. Huse, Phys. Rev. Lett. 70, 3147 (1993). 

^° J. Ye, S. Sachdev, and N. Read, Phys. Rev. Lett. 70, 4011 (1993). 

M.J. Rozenberg and D.R. Grempel, Phys. Rev. Lett. 81, 2550 (1998). 

D.R. Grempel and M.J. Rozenberg, Phys. Rev. B 60, 4702 (1999). 

J.V. Alvarez and F. Ritort, J. Phys. A 29, 7355 (1996). 
" L. Arrachea and M.J. Rozenberg, Phys. Rev. Lett. 86, 5172 (2001). 

P.W. Anderson and G. Yuval, Phys. Rev. Lett. 23, 89 (1969); G. Yuval and P.W. Anderson, Phys. Rev. B 1, 1522 (1970); 

P.W. Anderson, G. Yuval, and D.R. Hamann, ibtd. 1, 4464 (1970). 

Y.Y. Goldschmidt, Phys. Rev. B 41, 4858 (1990); T. Obuchi, H. Nishimori, and D. Sherrington, J. Phys. Soc. Jpn. 76, 
054002 (2007). 

L. De Cesare, K. Lukierska-Walasek, I. Rabuffo, and K. Walasek, J. Phys. A 29, 1605 (1996); T.M. Nieuwenhuizen and F. 
Ritort, Physica A 250, 8 (1998). 

T.R. Kirkpatrick and D. Thirumalai, Phys. Rev. Lett. 58, 2091 (1987); Phys. Rev. B 36, 5388 (1987); A. Crisanti, H. 

Horner, and H.-J. Sommers, Z. Phys. B: Condens. Matter 92, 257 (1993). 
1^ Y.L. Wang, S. Shtrikman, and H. Callen, Phys. Rev. 148, 419 (1966). 
2° R.B. Grifiiths, Phys. Rev. Lett. 23, 17 (1969); B.M. McCoy, ibtd. 23, 383 (1969). 

L. Arrachea and M.J. Rozenberg, Phys. Rev. B 65, 224430 (2002). 



